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We propose a non-equilibrium, hadronic kinetic model for describing the relative abundance of 
strange particles in ultra-relativistic heavy-ion collisions. The energy dependence of the multiplic¬ 
ity ratios of charged kaons and lambdas to pions are studied in detail. The pronounced peak in 
{K'^) /(tt’*') is conjectured to be due to the decreasing lifetime of a fireball with an increase in beam 
energy above 30 AGeV. Such a behaviour may be the consequence of a stronger stopping at lower 
beam energies. We study the dependence of the multiplicity ratios on energy and baryon density of 
the fireball and its total lifetime. 


I. STRANGENESS PRODUCTION IN 
NUCLEAR COLLISIONS 

The excitation functions of multiplicity ratios of kaons 
to pions are among the most intriguing results of the en¬ 
ergy scan at CERN’s Super-Proton-Synchrotron 0,11 
The sharp peak of the ratio (itr+)/(7r+) for collisions at 
beam energy of 30 AGeV—dubbed “the horn”—is barely 
reproduced by any theoretical calculation. The statisti¬ 
cal model expects a maximum in this region, which is 
due to the transition from baryon-dominated to meson- 
dominated matter Q. That maximum is, however, much 
broader than the observed sharp peak. Transport mod¬ 
els UrQMD and HSD ^ overpredict 7r+ production and 
predict a very modest bump which is lower than the mea¬ 
sured horn and lies at lower beam energy. The BUU 
model of shows no peaW structure at all. Three- 
fluid hydrodynamic model |g reproduces the K'^ and 
pion yields but overpredicts a multiplicity of K~ . Re¬ 
cently, a kinetic model which assumes a phase-transition 
to the deconfined phase has beeiiproposed and was suc¬ 
cessful! in reproducing the data Q- 

The data can be interpreted in the framework of the 
Statistical Model of the Early Stage (SMES) In 
this model it is assumed that the primordial particle 
production follows the prescription of statistical equilib¬ 
rium. The increase of the {K ^)/ratio at low colli¬ 
sion energies is a consequence of an enhanced strangeness 
content due to a larger energy deposit. The sharp de¬ 
crease of the ratio above the beam energy of 30 AGeV 
is connected with a first order phase transition between 
hadronic and deconfined phases. The decreasing part 
of the excitation function corresponds to the mixture of 
the two phases, while the flat dependence at high ener¬ 
gies is realized in the deconfined phase. In framework 
of SMES, this observation is taken as partial evidence 
for deconfinement in the nuclear collisions at threshold 
energy ^'snn « 7.6 GeV (Ebeam = 30AGeV). 

At this point we would like to pose two questions which 


provide the motivation for our study. Firstly, is the as¬ 
sumption of chemically equilibrated primordial produc¬ 
tion justified? This is not obvious at all. Analysis within 
a statistical hadronization model shows that the limited 
volume in proton-proton collisions causes suppression of 
strangeness once canonical statistics is used. On the 
other hand, a large volume allows for the use of grand- 
canonical statistics and leads to enhanced strangeness 
content But how could the incident nucleons be in¬ 
fluenced by volume effects in the very moment of their 
first interactions? Secondly, one cannot claim evidence 
of deconfinement unless all hadronic scenarios have been 
safely ruled out. We shall not touch here the question of 
early thermalization at the SPS. In this paper we inves¬ 
tigate if all hadronic (non-equilibrium) models are safely 
excluded after comparison with the data. 

We propose a model in which strangeness is pro¬ 
duced dynamically in reactions involving the constituent 
hadrons of the fireball. Two aspects are crucial for the 
final total amount of produced strangeness: the energy 
density and the time. The energy density is reflected in 
the temperature and an increase of temperature enhances 
the strangeness production rate. On the other hand, at 
least for a non-expanding system, if the initial amount 
of strangeness is below the equilibrium value, it will cer¬ 
tainly increase with time until it saturates. Thus a longer 
lifetime also means higher strangeness content. Based on 
these considerations we speculate that the initial increase 
of the (Ar+)/(7r+) excitation function is due to increasing 
energy density. Can the following sharp decrease of the 
excitation function be due to decreasing total lifetime of 
the fireball with increasing collision energy? 

The assumption of decreasing total lifetime is crucial. 
Sometimes, just the opposite is expected. Due to a higher 
initial energy density at higher beam energies it would 
take a longer time for the hreball to cool down to the 
critical density when it breaks up. Also, a larger mul¬ 
tiplicity at higher energies assumes a larger volume at 
freeze-out and thus requires longer lifetime. These argu- 
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ments, however, neglect nuclear stopping and the need to 
build up longitudinal expansion. For example, the proper 
time between impact and the freeze-out comes out more 
or less independent of the beam energy in the transport 
code HSD [2. Indeed, it seems reasonable to expect 
that stopping power or an ability for slowing down inci¬ 
dent baryons is stronger at lower energies than at higher 
energies. Then, a large part of the longitudinal expansion 
flow must be built up by the pressure at lower collision 
energies, and this takes time. At higher collision energies, 
fast expansion may be present from the very beginning 
and the hreball would reach a breakup stage earlier. We 
should stress that there is no direct measurement of stop¬ 
ping power as it results in rapidity distributions just after 
the primordial nucleon-nucleon collisions. Measured ra¬ 
pidity spectra reflect distributions at freeze-out, i.e., af¬ 
ter the possibly accelerated longitudinal expansion. We 
know that there is some stopping at RHIC since the net 
baryon spectra are peaked at a rapidity lower than ybeam 
im We also know that stopping is incomplete even down 
at SIS energies (few hundreds AMeV) as it was demon¬ 
strated by isospin tracing in collisions of 44 Ru and loZr 
M It is natural to expect that the relevant evolution 
scenario at AGS and SPS energies will be somewhere be¬ 
tween full stopping and re-expansion (Landau scenario) 
and boost-invariant non-accelerated longitudinal expan¬ 
sion (Bjorken expansion). 

In our model we shall adopt an ansatz for the time 
dependence of energy density and baryon density of 
the fireball. Our calculation is not based on any 
transport/cascade models or hydrodynamic simulation. 
Therefore we do not directly connect to details of micro¬ 
scopic structure of the matter or the equation of state. 
We gain, however, the freedom to construct as many 
different evolution scenarios as we wish, and we shall 
use this freedom to explore various scenarios, their im¬ 
pact on data, and identify those which are allowed by 
data. In particular, we shall focus on the dependence 
of total production of strange particles on the lifetime 
of the system. We present results for multiplicity ratios 
(Ar+)/( 7 r+), {K~)/{'k~), and (A)/( 7 r). It will turn out 
that a non-equilibrium hadronic scenario cannot be ex¬ 
cluded by this set of data. 


II. THE MODEL 

We shall only calculate ratios of yields and not the 
yields themselves. Therefore, it will be sufficient to know 
average densities of individual species if we assume that 
they come from the same volume. (Feed-down from reso¬ 
nance decays is also necessary, but we defer its discussion 
for later.) Applying later this approach to kaons (AT’*' 
and Ai°), we implicitly assume that kaons do not decou¬ 
ple from the system until the overall freeze-out. This is 
in line with a class of hydro-chemical models [Ullllll 
This assumption will be supported by the analysis of 
the kaon scattering rate and the mean free path in Sec¬ 


tion 

The evolution of kaon density in proper time is de¬ 
scribed by the master equation which can be derived from 


dpK d Nk ( 1 1 dNK 

The first term on the right-hand-side is due to the change 
of a fireball volume and includes the expansion rate. The 
expansion rate is the change of the density of any con¬ 
served charge due to expansion. For example the relation 


I dV _ I dpB 
V dr pb dr 

holds for the baryon number. The second term on the 
right-hand-side of eq. o stands for the change of a kaon 
number due to production or annihilation reactions. It 
can be divided into gain term and loss term 
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Terms in angular brackets are the cross-sections for reac¬ 
tions ij —*■ KX or Ki —> X which are multiplied by rela¬ 
tive velocities of incoming particles and averaged over dis¬ 
tribution of relative velocities. Here, i and j denote any 
single species with densities pi and pj , and X stands for 
any number of any species. The sums run, in principle, 
over all possible reaction channels. In practice they are 
truncated after the most significant contributions have 
been included. 

Using eq. the master equation for kaon density 
evolution can be written 

dpK ( 1 dV\ , , 

^gain - 7^1oss (4) 


In our study we shall evolve densities of Ar+, Ar*+, 
and K*° according to this equation. We shall explain 
how densities for other strange species are obtained in 
Section III Bl 


A. Expansion dynamics 

In order to calculate the production and annihilation 
rates we need the densities of interacting species and dis¬ 
tribution of relative velocities. These follow from the 
energy density £ and the densities of conserved charges, 
baryon number pB and the third component of isospin 
P 3 , which is also related to the electric charge. We 
shall assume thermal equilibrium and use £, ps and ps 
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to calculate temperature. The temperature determines 
the distribution of relative velocities. From the time- 
dependence of pb we shall obtain the expansion rate. 

We want to be able to explore a variety of possible 
expansion scenarios. We shall tune the amount of stop¬ 
ping which is directly connected with the time-derivative 
of the density at the beginning. Less stopping leads to 
larger negative time-derivative of ps because of the ini¬ 
tial stage expansion. Explored scenarios can also differ 
in the initial acceleration due to inner pressure. To in¬ 
clude these two effects in a simplest possible way, we use 
quadratic dependence as an Ansatz for the initial time- 
dependence of the energy density 

e(T) = £o(l — ar — br^) (5a) 

with £ 0 , a, and b being tunable parameters. In corre¬ 
spondence to this we use for the initial time dependence 
of conserved charge densities 

Piir) = pioil - ar - br"^)^ (5b) 


where the index i stands for any of the two conserved 
charges. The power 6 is used here as the simplest way 
of representing the equation of state. For example, in 
case of non-interacting freely streaming gas no work is 
performed by pressure and the energy density follows the 
same time-dependence as baryon density, 5 = 1. On the 
other hand, maximum pressure is present in a massless 
gas, in which ps oc 5 = 3/4. We shall write 


1 




1 + K 


o<k< 5 . 


( 6 ) 


The dependence of Hanbury-Brown/Twiss correlation 
radii on the transverse momentum of a pair indicates that 
at late times the fireball expansion follows a power-law 
in time 



(7a) 

'’•w = 

(7b) 


The power a is a model parameter. The value of a = 1 
corresponds to a one-dimensional boost-invariant expan¬ 
sion (Bjorken scenario). Due to transverse flow, a ac¬ 
quires values larger than 1. The parameter tq allows for 
a shift of the power-law prescription o to later times 
and leaves thus time for acceleration. 

Ansatz for initial times is matched with the pre¬ 
scription for late times o at some moment Tg. We re¬ 
quire continuous time dependence together with its first 
time derivative. In summary, the time dependence is 
given by 


£(t) = <1 

f £ 0(1 — OT — 6 t^) : 

T <Ts 


T > Tg 
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(8b) 


For the energy density we have parameters Eq, a, b, £g, 
Tq, Tg, a, and 5, out of which two are constrained by 
matching conditions at Tg. In addition to this, for both 
baryon and I 3 density there are pio and p'g which are 
connected by the requirement of continuity as well. 

The parameters which appear in eqs. m are not 
directly related to measurements and/or one may not 
have good feeling for what their reasonable values should 
be. Therefore, for practical application, we shall spec¬ 
ify the time evolution in terms of the following quantities: 
£0 initial energy density, 

tt total lifetime of the fireball until chemical freeze- 
out after which no change of number of strange 
hadrons occurs; this time must be bigger than Tg, 
Ef final energy density, 

Pb f final baryon density, 

P 3 f final density of the third component of isospin I 3 , 
Tq reference time for the power-law dependence in 
the late stage of evolution; in principle, it is the 
time (tt — tq) which would appear in expansion 
velocity gradients at freeze-out; (t — tq) is also 
the time measured by longitudinal HBT radius, 
a exponent of the late-stage power-law expansion; 
a = 1 corresponds to one-dimensional scaling ex¬ 
pansion, a = 3 means three-dimensional Hubble 
flow 

6 exponent characterizing the underlying equation 
of state, 

^max the maximum expansion rate, which is equal to 
the negative of maximum decrease rate of baryon 
density (—l/ps) (dps/dT) {TZmax is always pos¬ 
itive number); the maximum rate is achieved at 

Ts- 

We list in Appendix m relations between these physical 
parameters and those of parameterizations ©• 

The main parameters, which influence on the final 
strangeness yields we shall study below, are the initial 
energy density £0 and the total lifetime of the fireball tt- 
The final state parameters £/, pBf, P 3 f will be calculated 
from the thorough analyses of various particle ratios in 
the chemical freeze-out model ^3. Parameters tq and 
a will be chosen to agree with femtoscopy analyses na 
We shall choose a value for 5 and fix it and i?niax can be 
varied only in very limited range and has little influence 
on the results. 

The expansion rate needed in equation 0 is obtained 
by taking time-derivative of eq. Il8bll . Using eq. 0 we 
have 


1 ^ 


(5(a+2&r) 
l — aT — br-^ 


r <Ts 

T>Ts 


(9) 


For illustration we plot in Fig. Q] the time dependence 
of the energy density and baryon density and the temper¬ 
ature. The curves correspond to the scenario which re¬ 
produces experimental data on (iF+)/(7r+), {K~)/{tt~), 
and {A)/{tt) ratios for 30 AGeV collisions as it is shown 
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FIG. 1: Left panel: the energy density (left scale) and the 
baryon density (right scale) as a function of time for a scenario 
that reproduces the data from 30 j4GeV collisions. Right 
panel: the corresponding time dependence of the tempera¬ 
ture. 


end. This is driven by the requirement of a rather low ini¬ 
tial energy density. Since the final density is fixed by the 
data, a longer power law tail leads readily to very high 
initial energy densities. We emphasize that the proposed 
evolution does not contradict to observations related to 
the freeze-out stage of a fireball expansion reconstructed 
via single-particle spectra and femtoscopy. 


B. Chemical composition and reactions 


For the sake of averaging over relative velocities, we 
shall assume that the momenta are distributed according 
to Boltzmann distribution 


in Fig. El with parameters listed in Tables ^ and cni 
Time evolutions in other scenarios are qualitatively very 
similar to the one presented here. It turns out that the 
shape of the curves is almost specified by the initial and 
final energy density and the lifetime. In fact, most of the 
time the fireball spends in the regime of accelerated ex¬ 
pansion while the power law is realized only towards the 


mip) oc exp . (10) 

Thus we make an assumption of thermal equilibrium, 
though chemically the system will be treated as non- 
equilibrated. The averaged cross section is then obtained 
as 


Vzjcrf^) = 


I^dxcrf^ (x) Ki (f) [a:^ - (m^ -b rrij)^] [a;^ - (m, - 


4mf m]TK2{m,/T) 


( 11 ) 


where Ki% are the modified Bessel functions and = 
max(mi -I- m^ , X)finai”^a) reaction threshold. 

The following reaction channels producing kaons were 
taken into account 
-kN ^ KY ttN 

ttA KY ttA 

ttY ^ KE NN 

TTTT KK NA 

np <—!■ KK NA 

np <—!■ KK* AA 

pp ^ KK NN 

K* ^ Ktt 

In order to keep the detailed balance we also included 
the inverse reactions for all channels with two particles 
in final state. Those processes with three and more final 
state particles are rather suppressed due to high thresh¬ 
olds and smaller phase space. By the same reason we ne¬ 
glect strange antibaryons. The error we thus introduce is 
small, because the matter is rather baryon-dominated in 
the investigated energy domain. The ratio A/A is about 
10% at the highest SPS energy, so our discrepancy will be 
at most of this order. This is acceptable for a schematic 
model like the present one. 

Only production of K = {K'^,K0) and K* = 

(A'*+,iL*°) mesons will be calculated from master equa¬ 


tion 0 by including the above reaction channels. All the 
explicitly calculated species contain strange antiquark. 
They are produced being accompanied by strange and 
multistrange baryons, and K and K* mesons contain¬ 
ing a strange quark. In principle, we could trace density 
evolutions of those species kinetically. In practice, how¬ 
ever, the reactions which do not create strange quarks 
but rather rearrange them in different hadrons are very 
quick. Therefore with a good approximation we can as¬ 
sume that all species containing strange quarks are in rel¬ 
ative chemical equilibrium. In our calculations we include 
K-, K°, K*-, K*°, A, E, A(1405), A(1520), S(1385), 
S, and n. 

For the non-strange sector we assume that all species 
are chemically equilibrated. We included into calcula¬ 
tions all non-strange mesons with masses below 1.5 GeV 
and baryons up to 2 GeV. 

Since we follow evolution of every isospin state sepa¬ 
rately, we do not use isospin-averaged parameterizations 
for the cross sections. Appropriate cross section for every 
isospin channel is used instead. We list the used param¬ 
eterizations in Appendix IBI 

For the calculation of densities of individual species 
and the evaluation of the average in eq. m we thus 
need the temperature T, and chemical potentials pB and 


NKK 

NA - 

-4 NNKK 

NKK 

AA - 

-4 NNKK 

KNY 

NYK 

AKY 

AYK 

AKY 

NN - 

NNKK 


( 12 ) 
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/i 3 connected with the baryon number and the third com¬ 
ponent of isospin, respectively. We also need phase-space 
occupation factors for kaons 7 if+, 7 /<-o, 7 _ff+, and 
as well as one common factor for species with S < 0, 75 . 
On the other hand, the time evolution is formulated in 
terms of e, ps, pa, and kinetically calculated px+i Pk°, 
Pk*+, and pk*o. We thus need to write down relations 
connecting these two sets of quantities. First, we note 
that the overall strangeness neutrality dictates that the 
density of strange quarks 

PS = \Si\pi = px+ + Pk° + Pk*++ ( 13 ) 

i, S<0 


where the sum should run over all species containing 
strange quark, though practically we include only those 
mentioned above. Then, in Boltzmann approximation, 
we have the following relations between thermodynamic 
quantities and densities 


e = 


27r^ 


E 9 ^\rnlT^ 


(14a) 


i.S=0 




27r^ 


E \ ISil 2 rr2 
gz Xz 7s 


i,S<0 


Pb = 


+ ^ E giXzlzmjT'^ 

kaons 

^ (14b) 


i,S=0 


+ ^ E ^*53A^7^®‘'mfTX2 (^) 


2,5<0 


P3 = 


^ Y (14c) 

^ 2 ,S —0 


2,S<0 


+ i E h,zgzXzlzrnjTK2{^) 


2 , kaons 


2,S<0 


^ 2^ "^^2 (^) (14d) 


and for all kaon species 
1 

27r^ 


1 (h,zP3 \ 2rj.r^ fmi\ 

Pz = iTZ2 1 ) 7* w, T K 2 — J , 


T 


i = K+, K°, iW+, K*° . (14e) 


These relations can be inverted numerically and T, ps, 
P 3 , 7 S) and the four y^’s can be determined. 

In equations 03, Pi is the spin degeneracy of species 
i (different isospin states are accounted for separately). 
The fugacity 


Ai = exp((BjpB-|-/3,ip3)/r) , (15) 

Bi, I^^i, and Si are baryon number, third component of 
isospin, and strangeness, respectively. 


C. Kaons as a part of the system 


Kaons have rather small cross sections for interactions 
in baryon-rich matter and therefore their scattering rate 
(mean number of collisions per unit of time) is not large 
and the mean free path is long. In our simulations we 
calculate the scattering rate in a similar way as the an¬ 
nihilation rate in eq. 


7?.scatt — 'y 

i 


Pz 

1 + Sxi 


(16) 


In determining 7?.scatt we take all processes that are in¬ 
cluded in kaon annihilation and add the elastic scattering 
off protons and neutrons with the cross sections given by 
eqs. (IB44II and IIB45II . 

Kaon scattering rates for three representative expan¬ 
sion scenarios are plotted in Fig.|21 They are typically of 
the order 0.1 fm“^; somewhat larger at higher density in 
the beginning and smaller at later times. With thermal 
kaon velocity of about 0.7 c this leads to a mean free path 
of 7 fm. Even if the kaon moved with light velocity, like 
some primordially produced kaons, the mean free path 
would still come to about 10 fm. This length is compa¬ 
rable to the size of a fireball. Would the kaon mean free 
path be much bigger than the size of a fireball, one could 
assume kaons to escape from the system right after the 
production event without any further scattering. Oppo¬ 
sitely, a mean free path much smaller than the fireball 
size would indicate that kaons rescatter intensively and 
stay in thermal equilibrium with fireball medium. Real¬ 
ity is somewhere in between these two limits. In any case, 
they do not simply escape from the fireball. In our type 
of model we can only make two extreme assumptions: de- 
conple the produced kaons from the system completely or 
treat them all as a part of the system. Due to a chance of 
absorption the latter appears closer to reality. For prac¬ 
tical reason, in calculation of kaon scattering and annihi¬ 
lation rates thermally equilibrated velocity distributions 
have been assumed. 

Let us also note that the inverse slope parameters T* 
of kaon spectra, together with pions and protons seem to 
follow the empirical prescription T* = T -\-m{vtj\ where 
(ut)^ is average transverse expansion velocity |l9l| . This 
suggests that all these species freeze-out kinematically at 
the same time. 
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FIG. 2: Comparison of the kaon scattering rate calculated according to Eq. m and the expansion rate introduced in Eq. 0 
as a function of time for three scenarios which fit the data at: 11.6 ^IGeV (left), 30 AGeV (middle), and 158 AGeV (right). 


The influence of our simplification is twofold. Firstly, 
we may overestimate the amount of kaons that are an¬ 
nihilated. Secondly, the energy which the kaons carry is 
not taken out of the system as it would be so if kaons 
decouple. 

The overestimation of kaon annihilation would intro¬ 
duce only small error if the annihilation rate is small 
compared to the production rate. The rates are plotted 
for three typical scenarios in Fig. |21 (same scenarios as 
in Fig. n. The annihilation rate grows since the kaon 
abundance increases. Originally it is smaller than the 
creation rate by an order of magnitude, at late times by 
a factor of 2. From the fact that up to 50% of final state 
direct kaons are produced primordially and not due to 
kaon creation modelled here (see Fig. CH), we conclude 
that a decrease of the (small) annihilation rate would not 
change our results too much. At most, it would slightly 
shorten the lifetimes necessary to describe the data. 

In order to estimate a possible error in the total en¬ 
ergy of the fireball we plot in Fig. 0 the ratio of the 
energy carried by thermal kaons to the total energy of 
the system. At the beginning, at most 6 % of energy is 
contained in kaons. This is the part of energy we should 
take away if we treat primordial kaons as non-interacting. 
Since energy density goes roughly with the fourth power 
of temperature, T would not change much and the pro¬ 
duction rates would stay practically unchanged. In gen¬ 
eral, subtracting a few per cent from the energy density 
would have small influence on chemical composition be¬ 
cause at the same time we would not include kaons and 
the available energy would be distributed among smaller 
number of degrees of freedom, thus changing the temper¬ 
ature only little. 

It is possible to refine our treatment of kaons as a part 
of the fireball by comparing the kaon scattering rate to 
the expansion rate of the fireball, cf. Fig.|21 At the be¬ 
ginning of fireball evolution, expansion rate is lower than 
the scattering rate, thus we say that kaons are thermally 
equilibrated. When the expansion rate becomes compa¬ 
rable or bigger than the scattering rate, this means that 
the density of fireball drops considerably before a kaon 
had a chance to scatter. Then kaons are assumed to de¬ 


couple from thermal equilibrium. Decoupling is a gradual 
process but we cannot treat it so in the framework of our 
model. We say that the kaons decouple at once when 

^ >07^,catt. (17) 

Pb dr 

When this condition becomes fulfilled, we hx the tem¬ 
perature for kaons in eq. 114e|l . We thus assume that the 
kaons keep a higher temperature than other species, be¬ 
cause they decoupled from the thermal bath earlier. By 
default, the constant 9 is set to F though no exact value 
can be specified by theory . Variations of 9 turn 

out to have negligible influence on the results. 

In Fig. [21 we demonstrate that the decoupling accord¬ 
ing to this prescription happens only towards the end of 
expansion. 

D. A summary of the evolution procedure 

For an easier overview, we summarize the algorithm 
for time evolution. Suppose, at time r we know all quan¬ 
tities: e, Pb, P3, all kaon densities, and thermodynamic 
quantities T, pB, P 3 , 7S, and all y^’s. We proceed to 
time r -|- dr in following steps: 

1. Calculate new densities of , K^, A'*+, K*^ from 

eq. 0). 

2. Calculate total kaon scattering rate TT-scatt from 

eq. (tTBl) . 

3. Obtain e, ps, ps from prescription ® . Also obtain 
expansion rate from eq. ®. 

4. Determine ps from eq. ca. 

5. Decide by the use of inequality 03 if kaons are 
thermalized according. 

6 . Obtain T, ps, ^ 3 , 7s, and y^’s from numerically in¬ 
verting relations 03- If kaons are not thermalized, 
fix their temperature in eq. 114ell . 
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FIG. 3: Comparison of total kaon creation and annihilation rates as a function of time for three scenarios which fit the data 
at: 11.6 ^GeV (left), 30 AGeV (middle), and 158 AGeV (right). 





FIG. 4: Ratios of the energy contained in kaons to the total energy of the system as functions of time for three scenarios which 
fit the data at: 11.6 AGeV (left), 30 AGeV (middle), and 158 AGeV (right). 


7. Calculate the density of any desired species via 

= Ai mj T K2 (18) 

for non-strange species, and 

= K 7^' mjTK2 (^) (19) 

for species with 5” < 0 . 

8 . Continue to next time step by going to step 1. 

In the outlined setup of the model there are differ¬ 
ent handles to tune the final density and those of 
K~ and A. All these, normalized by pion densities, will 
be compared to data. Since the production of is 
calculated explicitly, its density depends on temperature 
and time. If we can reproduce the {K^) /( 7 r+) ratio, this 
means that we have gotten the total strangeness pro¬ 
duction right. The strange quarks are then distributed 
among AT”, A, and other species according to the tem¬ 
perature and chemical potentials. Thus the key to si¬ 
multaneous fit to {K~)/{-k~) and (A)/( 7 r) ratios is the 
correct value of the temperature. 


E. Final state and feed-down from resonance 
decays 


Chemical composition of the final state in Au-|-Au col¬ 
lisions at beam energy 11.6 AGeV and 

Pb-fPb collisions at beam energies 30 AGeV 40, 80, 
and 158 AGeV Q have been analyzed in the framework 
of the statistical hadronization model complemented by 
strangeness suppression factor 75 by Becattini and 
collaborators [ig. We use their results on temperature, 
chemical potentials and 75 in order to characterize the 
state of chemical freeze-out which we aim for. Thus ob¬ 
tain the final state value of the energy density, B and 
/a densities (Table D). These are three of the parameters 
which specify the time evolution. Note that even fix¬ 
ing the final state densities does not guarantee that we 
arrive at the correct final chemical composition because 
final temperature depends on the amount of produced 
strangeness which is calculated kinetically. On the other 
hand, it is important to realize that if we end up in a 
state with the correct amount of kaons and the correct 
final energy and baryon densities guaranteed by construc¬ 
tion, then we also reproduce the temperature as inferred 
in chemical freeze-out fits in |16| and therefore all ratios 
of abundances, i.e. also those involving only non-strange 










TABLE I: Thermodynamic quantities and densities charac¬ 
terising the chemical freeze-out state in collisions at different 
energies. Values of T, ps, pa, fj,s, and 7 s from serve as 
an input in reconstructing the final state values of e/, ps f 
and p 3 f. 


Abeam [AGeV] 

11.6 

30 

40 

80 

158 

T [MeV] 

118.1 

139.0 

147.6 

153.7 

157.8 

PB [MeV] 

549.1 

423.1 

375.1 

293.6 

243.7 

P3 [MeV] 

- 11.8 

- 11.0 

-10.3 

- 8.2 

-7.1 

PS [MeV] 

117.5 

99.1 

90.5 

69.7 

59.2 

7s 

0.652 

0.938 

0.757 

0.73 

0.843 

£/ [GeV/fm®] 

0.132 

0.173 

0.203 

0.194 

0.198 

PBf [fm-®] 

0.086 

0.087 

0.091 

0.068 

0.058 

P3f [fm"®] 

-0.0083 

-0.0097 

- 0.0102 

-0.0079 

-0.0070 


TABLE II: Mean multiplicity (K^) NN from nucleon-nucleon 
collisions and ratios {K'^)NN/{h~)NN extrapolated to a mix¬ 
ture of protons and neutrons corresponding to Au-|-Au or 
Pb-|-Pb collisions at indicated beam energies. 



{K^)nn 

(h~)NN 

Au+Au @11.6 AGeV 

0.0607 ± 0.0252 

0.0581 

Pb-tPb @ 30 AGeV 

0.143 ±0.027 

0.0821 

Pb-tPb @ 40 AGeV 

0.163 ±0.036 

0.0842 

Pb-tPb @ 80 AGeV 

0.218 ±0.059 

0.0880 

Pb-bPb @ 158 AGeV 

0.272 ± 0.056 

0.0883 


species. If we then compare calculated density of one 
of the species to the total measured multiplicity of that 
species we could infer the volume and thus multiplicities 
of all other species. 

A significant portion of final state pions stems from 
decays of resonances. These must be included when ob¬ 
taining ratios of multiplicities. This is done by deter¬ 
mining the number of resonances and multiplying by the 
average number of pions produced in a decay of one reso¬ 
nance. Also, there is feed-down to kaon production from 
K* decays and to the number of A’s from other strange 
baryons. A list of all contributions is provided in Ap¬ 
pendix El 


F. Initial state 

Now we have to specify the initial conditions for the 
evolution equation 0 . Note that strangeness is also pro¬ 
duced in primordial collisions of incident nucleons. The 
amount will be extrapolated from strangeness abundance 
in nucleon-nucleon interactions. We shall equate the ra¬ 
tio of K'^ density and the density of negative hadrons 
to the ratios of multiplicities {K'^)jsin / {h ~)nn extrap¬ 
olated from nucleon-nucleon collisions. The values are 
summarized in Table m Density of is then deter¬ 
mined as pko = Pk+ exp(—/X3/T). No K*'s are assumed 


in the initial state. Species with S < 0 must balance the 
total strangeness to 0 and are set into relative chemical 
equilibrium. We describe the choice of the initial state in 
detail in Appendix [0 


III. RESULTS AND DISCUSSION 

We calculated the resulting (Ar+)/(7r+), {K~) / {'k~) 
and (A) /(tt) ratios for many different evolution scenarios 
for the beam energies 11.6 AGeV (Au-l-Au), and 30, 40, 
80, 158 AGeV (Pb-I-Pb). Thus we explore the region of 
the peak in (Ar+)/(7r+). 

At every beam energy the investigated scenarios differ 
by the initial energy density and total lifetime of the fire¬ 
ball. The final co nfig urations listed in Tableware fixed 
from the analysis [2l ■ The range of initial energy densi¬ 
ties goes up to 3.0 GeV/fm^. This seems to surpass the 
critical value for deconfinement indicated by lattice cal¬ 
culations. We have two comments here, however: first, 
the system under investigation is out of equilibrium and 
it is a question whether results of equilibrium statistical 
physics can be applied here. The relaxation time is not 
known. Second, an uncertainty in determination of the 
critical temperature was reported at the Quark Matter 
2005 conference leading to a value higher than one 
previously used. 

The parameter tq is chosen such that it is smaller than 
the total lifetime by 7 fm/c, tq = tt — 7fm/c. Such a 
value is motivated by measurement of the longitudinal 
HBT radius 77;, which—in a Bjorken scenario—indicates 
that the thermal freeze-out happens about 8-9 fm/c after 
the reference tq. As we deal in our study with chemical 
freeze-out, we set tt — tq a, little smaller. 

Since we know from transverse momentum spectra and 
HBT radius anlyses that there is considerable transverse 
expansion, we fix a = 2. The equation of state is set by 
a modest value ot k = 1/6, or J = 6/7. 

The last parameter to specify is 77.niax- If we fix the ini¬ 
tial and final densities and the total lifetime in which the 
fireball evolves between these two states, there is only a 
limited range of values the maximum expansion rate can 
assume. For every pair of chosen initial energy density 
and lifetime we explored the two limiting values of TT-max- 
It turned out that the differences due to selection of TT-max 
in the studied ratios were of the order 10%. 

The results are plotted in Figs. I5I9I The shown 
results for energies 30 AGeV and upward were obtained 
for the highest possible TT-max- If scenarios with slowest 
possible TT-max are used, the lifetime needed to reach the 
same result is prolonged by 1-2 fm/c. For the lowest 
studied energy we used the lowest TT-max scenarios. 

In Figs. l5l9l we clearly observe that the resulting multi¬ 
plicity ratios depend crucially on the total lifetime of the 
fireball. Dependence on the initial energy density is less 
important. The latter is more strongly pronounced at 
lower energies, where the baryon density is higher. Since 
the dependence on the energy density is so weak, we do 
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FIG. 5: Ratios {K^)/{K~)/{'k~), and {A)/{'k) as a 
function of the total lifetime of the system, calculated for 
Au+Au collisions at beam energy of 11.6 AGeV. Different 
curves correspond to different initial energy densities. Hori¬ 
zontal lines indicate la intervals around the measured data. 


FIG. 7: Ratios {A'^)/(7r^), {K~)/{'k~), and (A)/{7r) as 
a function of the total lifetime of the system, calculated 
for Pb-|-Pb collisions at beam energy of 40 AGeV. Different 
curves correspond to different initial energy densities. Hori¬ 
zontal lines indicate la intervals around the measured data. 
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FIG. 6: Ratios (JG^)/(7r^), {K~)/{'k~), and (A)/{7r) as 
a function of the total lifetime of the system, calculated 
for Pb-pPb collisions at beam energy of 30 AGeV. Different 
curves correspond to different initial energy densities. Hori¬ 
zontal lines indicate Icr intervals around the measured data. 


FIG. 8: Ratios {K^)/{'k^), {K~)/{'k~), and {A)/{'k) as 
a function of the total lifetime of the system, calculated 
for Pb-pPb collisions at beam energy of 80 AGeV. Different 
curves correspond to different initial energy densities. Hori¬ 
zontal lines indicate Icr intervals around the measured data. 


not expect a major change in the results if parameter¬ 
ization for time dependence of the energy density was 
changed. 

In order to quantify the agreement between our calcu¬ 
lations and experimental data, for every examined sce¬ 
nario we calculate the standard measure 

^2 ^ (theory, - dataQ^ 

^ error? ’ 

i ^ 

where the sum runs over the three measured data points 
for every beam energy: (ifr+)/(7r+), {K~) j{A)/{- k). 


The resulting values are summarized in Fig. irni For 
beam energies above 30 AGeV the best agreement with 
data is indeed obtained if the lifetime decreases as a func¬ 
tion of beam energy. At the highest AGS energy, our 
model fits data best if the lifetime is shorter than at 
30 AGeV. In some sense we have translated the non¬ 
trivial excitation function of {K^)/{tt^) into the depen¬ 
dence of the lifetime on the collision energy. A maxi¬ 
mum of the lifetime at 30 AGeV could be an indicator of 
some change in evolution dynamics, e.g. a soft point in 
the equation of state leading to a low pressure gradient 
and a prolonged lifetime. However, we can not insist on 
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[fm/c] 


FIG. 9: Ratios {K^)/{K~)/ {'k~), and (A)/(7r) as a 
function of the total lifetime of the system, calculated for 
Pb+Pb collisions at beam energy of 158 AGeV. Different 
curves correspond to different initial energy densities. Hori¬ 
zontal lines indicate la intervals around the measured data. 
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[fm/c] 

FIG. 10: The quality-of-fit measure , as defined in eq. 
as a function of the total lifetime. Different panels correspond 
to different beam energies and various curves to different ini¬ 
tial energy densities. 


TABLE III: Parameters of the simulations which lead to re¬ 
sults compared to data in Fig. mi In the lower two rows we 
compare the chemical freeze-out temperature T obtained by 
Becattini et al. IT^ with the teperature Tf we obtain in our 
simulations. 


Fbeam [AGeV] 

11.6 

30 

40 

80 

158 

£0 [GeV/fm®] 

1 

1.5 

2 

2.25 

2.75 

TT 

[fm/c] 

25 

25 

20 

15 

15 

TZmax [fm 

0.286 

0.305 

0.333 

0.379 

0.374 

T 

[MeV] 

118.1 

139.0 

147.6 

153.7 

157.8 

Tf 

[MeV] 

114.7 

134.1 

143.3 

149.3 

153.6 


such conclusion here, since comparison with data cannot 
clearly exclude that a lifetime at beam energy 11.6 AGeV 
is as long as at 30 ^GeV. A slight change in the depen¬ 
dence of tt on the collision energy, could be expected if 
other time evolution parameterizations is used. 

In Table Cni we present the most favorable sets of the 
evolution parameters. We choose those sets where the 
fireball lifetime does not increase with an increasing beam 
energy. Gorresponding final densities are given in Tabled 
Our results for the ratios (Ar+)/(7r+), {K~)/{tt~), and 
{A)/{tt) are shown in Fig. In Table HTTI we see that 

the final state temperature agrees rather well with the 
results of chemical freeze-out fits by Becattini et al. . 
Note that we did not perform a calculation at the beam 
energy 20 AGeV due to lack of final state analysis from 
this collision energy in . 

Table mil and Fig.^Jare the main results of our study. 
In the subsequent subsections we illustrate the time evo¬ 
lution of various densities and rates. 


A. The evolution 

Evolution of kaon density in Pb-l-Pb collisions at 
30 AGeV is illustrated in Fig.El(left panel). The growth 
of pk+ is soon overrun by the decrease due to expansion. 
To cancel the expansion effect we have to normalize the 
kaon density by the baryon density. The latter, being 
density of a conserved charge, changes due to expansion 
only. 

The ratio Pk+I Pb (middle panel in Fig. IT^ shows a 
steady increase in time. The ratio of pii+ / Pt^+ (right 
panel in Fig. PHl finishes with rather high value which 
does not correspond to the measured (Ar+)/(7r+). The 
final result which does reproduce the data includes, how¬ 
ever, feed-down from resonance decays which contribute 
largely to pion production. 

In Fig. ll 3l we show the density evolution for some other 
species. When the densities are normalized to baryon 
density in order to get rid of the effect of expansion we 
clearly see that due to cooling the relative abundance of 
pions and nucleons as the lightest mesons and baryons 
increases with respect to other species. Normalized K~ 
and A densities rise continuously in time. 
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FIG. 11: Comparison of our results with data uu Calculations were performed with parameters listed in Tables ||] and ITTTI 



FIG. 12: Left panel: evolution of the kaon densities. Middle panel: evolution of the kaon density divided by the baryon number 
density. Right panel: evolution of the ratio Pk+Ip^+. 



FIG. 13: Densities of p, p"*", A, and K~ as functions 
of time for the scenario which reproduces 30 AGeV data. In 
the lower panel the densities are normalized to the baryon 
density. 


B. The rates 


Kaons are produced in various reactions. In Fig. lldi we 
show the production rates of different reactions as func¬ 
tions of time. The three shown examples refer to scenar¬ 
ios which reproduce data for 11.6, 30, and 158 AGeV in 
Fig. m As expected, due to an increase of the energy 
density the production rates grow with increasing beam 
energy. In all cases, the dominant contribution is due 
to TTTT reactions. A very important contribution comes 
from reactions of pions with hyperons. Recall that the 
corresponding cross section is not known experimentally 
and we chose a constant matrix element for this reaction. 
This introduces some uncertainty into our quantitative 
results. Note however that since the contribution from 

































12 



FIG. 14: Production rates of due to selected processes as a function of time for three scenarios which fit the data at: 
11.6 ^GeV (left), 30 AGeV (middle), and 158 AGeV. 


Try reactions has almost always the same relative impor¬ 
tance in comparison to tttt reactions, we do not expect 
any qualitative change of our results under modification 
of ttY cross sections. 

There is a clear change when moving from lower, 
baryon-dominated, energies toward higher, meson- 
dominated energies. While at 11.6 AGeV the second 
largest contribution comes from ■kN reactions, the con¬ 
tribution from this type of reactions decreases with in¬ 
creasing beam energy and actually becomes even lower 
than Tip contribution at 158 AGeV. 

We study the meson-meson production rates together 
with TriV, ttV and iVA rate as a function of temperature 
in Fig. cni As mentioned, the difference between the 
three different panels in in baryon density. The tttt rate 
becomes relatively more and more important. As it is 
outlined in Annendices IH 7l a,nd lR 81 we took a great care 
in determining the irp and pp cross sections. Note that we 
do not treat all meson-meson cross section as being the 
same, in contrast to 14,511 . The resulting suppression in 
TTp channel is due to p-wave suppression of the reaction 
Tip ^ (j) ^ KK and a higher threshold in the chan¬ 
nel Tip KK*. The later cross section was determined 
from the scattering amplitudes following from the covari¬ 
ant solution of the Bethe-Salpeter equation with the 
leading chiral order Weinberg-Tomozawa term as an in¬ 
teraction kernel. The untarization effects are found to be 
important leading to a smaller cross section than what 
would be inferred from tttt cross section. The pp rate is 
suppressed thermally due to large mass of the p mesons, 
and the resulting contribution is very low. 


IV. CONCLUSIONS 

In the introduction we motivated our study by the 
question whether it is possible to rule out any hadronic 
interpretation of the observed excitation function of 
(Ar+)/(7r+), {K~)/{tt~), and (A)/(7r) ratios. We pro¬ 
posed a non-equilibrium hadronic model which is able to 
reproduce the data by virtue of varying the total lifespan 


of a fireball at different energies. 

Thus one has to do more in trying to exclude this kind 
of models. Another indications of non-trivial effects in 
the region of the {K~^) /(7r+) peak is a step in kaon spec¬ 
tra inverse slopes, and a change in the pion yield per 
participant. The model must be tested on these observ¬ 
ables, as well. 

Our findings are based on the assumption that the life¬ 
time of a fireball becomes shorter as one moves to higher 
beam energies. We want to stress that most of the infor¬ 
mation we have about the space-time evolution of a fire¬ 
ball comes from hadronic spectra and correlations, which 
are formed at the freeze-out stage of fireball evolution. 
This is just the final state to which many possible evolu¬ 
tion scenarios may lead. One of the future tasks should 
be a careful check how hadronic observables calculated 
within our model fit the data. This is rather involved 
project. Although we were lead by the data in formu¬ 
lating our model, we refrained from such detailed and 
careful fits here. 

The whole duration of fireball evolution is reflected in 
penetrating probes like dilepton spectra. Their measure¬ 
ment was improved dramatically in the last years I3EJ 
113 ■ It seems to be crucial, therefore, to check our model 
assumption about the lifetime against the dilepton spec¬ 
tra. 

Our approach addressed here only a part of the whole 
problem by asking how do different fireball evolution sce¬ 
narios influence the data. We ignored the question how 
a specific evolution scenario follows from the microscopic 
properties of the created matter? This question must be, 
of course, addressed in the future if the model survives 
all experimental tests. On the other hand, if we succeed 
to exclude the model on a set of data, there will be no 
need to answer this question. This approach might be 
simpler to begin with. 
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T [GeV] T [GeV] T [GeV] 

FIG. 15: Production rates of from selected reaction processes as a function of temperature for three scenarios which fit 
the data at: 11.6 ^GeV (left), 30 ^GeV (middle), and 158 AGeV. 
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APPENDIX A: PARAMETERS OF THE TIME 
DEPENDENCE ANSATZ 


Here we express the parameters of eqs. © in terms 
of eo, tt, £/, PBf, Pa/, ^ 0 , a, and TT-max, which were 
introduced in Section III Al 

The easiest to obtain are £q, p^q, p'^q 

4 = £/ (.TT - ro)“'^^ (Ala) 

p'bo = Pb/(7-t-to)“ (Alb) 

Pso = P3/(tt-to)“ (Ale) 

Since the expansion rate at time Ts must be T^max, from 
eq. © we see 




a 


n 


max 


+ ro. 


(A2) 


TABLE IV: Parameters of the time dependence 0 corre¬ 
sponding to the energy density, hreball lifetime and 77max 
given in Table irm 


E-beam [AGeV] 
a X 10"^ 
b X 10'‘ [fm-2] 


Ts [fm/c] 


11.6 30 40 80 158 

142.6 5.482 4.963 10.77 8.448 • 10 

7.752 13.29 21.55 41.60 44.70 

25.00 24.56 19.01 13.28 12.98 


-2 


which lead to 


a 


b 


2 


Ts 

_££ ( ‘^ \ ^ ^ f rr-To 

£o \ts S{ts - To) J \Ts- To 

£/ / J_ o \ / Tt-Tq 

So \T^ Sts(ts - To)) \ts- To 
1 


ct/S 

(A4a) 


Oi/S 


(A4b) 


For B and I 3 densities we obtain from matching at Ts 


Pbq 

P30 


PBf f Tt-Tq Y^ 

(1 - aTs - br^)^ \Ts-To ) 

_Pa/_ /D—“ 

(1 - UTs - br^)^ \Ts - To ) 


(A5a) 

(A5b) 


In Table rm we present the values of parameters a, b 
and Ts corresponding to the energy density, fireball life¬ 
time and TT-niax given in Table ITTTI 


Parameters n and b are obtained from matching at time A.PPEN^DIX^ B* CROSS SECTIOI^S FOR ICA.ON^ 
Ts PRODUCTION 


£ 0(1 - UTs - brf) = 
£ 0(1 -f 2 bTs) = 


(ts - To)°‘/^ 
S{ts - 


(A3a) 

(A3b) 


In this section we specify the cross sections of all re¬ 
actions with strangeness production which are relevant 
at the collision energies under consideration. We disre¬ 
gard a possible modification of cross sections in dense 
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FIG. 16: Total cross section of the ^ K^A, n'^p —> and tv p ^ K'^'E reactions: solid line calculated with 

eq. IIB3II . dots are the experimental data from Ref. liil^ 


and hot nuclear matter and accept their vacuum forms. 
This approximation complies with our assumption above 
that properties of all particles do not change in medium 
either. 

For processes with two particles in the final state the 
inverse reactions are included in the calculation, as well. 
The cross section for the inverse reactions follows from 
phase-space considerations as 


where Ji and rrn are spins and masses of the participating 
species, and Pcm is the center-of-mass momentum defined 
as 


Plmis,ma,mb) = 


[s — (mf -I- m^)]^ — 4 m\ 


4s 


(B2) 


< 734^12 (-s/s) 


_ ( 2 J 3 -I- 1)(2J4 -I- 1) pl^{s,mi,m2) 
(2Ji -I- 1)(2J2 + 1) Pcm('S, W3,TO4) 

X 0'i2^34(-\/s) , (Bl) 


In all parameterizations below, energy units will be 
GeV and cross sections are measured in millibarns. 


I 

1. Reactions of tt -|- 


Cross sections for ttN K'^A are all related by isospin symmetry M to 7r°n ^ K^A, which we take from 

2a( A - K+A) = a{Tv+n ^ K+A) = ^ K^A) = 2a(Tr% - iF°A) = ■ (B3) 

In Fig. 1161 we compare this parameterization of the 7r'''p —> K^A reaction with the available experimental data. 
,l,From we get the ttN cross sections 


a(7r~^p —> = a{'K n —*■ ) 

cr( TT°p K+YP ) = cr( 7 r°n ^ K° Y° ) 
a{'K~p K^Y~) = cr(7r+n ^ K^Y'^) 


0.03591(Vs- 0.1594(Vi - 1.688 )°°i°'^® 

(75- 1.890)2-^0.01548 ”^ (Vs - 3.000)2-h 0.9412 ’ ^ > 

0.003978(7s- 1.688)°-5®4® 0.04709(7 s - 1.688)2-1650 

(7s - 1.740)2 -h 0.00667 (7s - 1.905)2 -h 0.006358 ’ ^ ’ 

0.009803(7i- 1.688)°®°2 i 0.006521(7i - 1.688)1-4728 

(7s- 1.742)2 -1-0.006583 (7s - 1.940)2 -i-0.006248 ’ ^ ’ 


a{TT'^n —*■ K'^Y'^) = a{TT p —> iF*'!]!’) 


a{TT^n —*■ K^Y ) = a{Tv^p —*■ K'^Y^) 


0.05014(7i- 1.688)1-2878 
(7s- 1.730)2 -h 0.006455 ' 


(B7) 


The quality of these parameterizations is demonstrated in was found that under assumption that main contribution 
Fig. Elfor TT+p ^ K'^Y~^, and 7r“p —s- K'^Y~ reactions. 

Reactions ttN ^ NKK were studied in Ref. [s^. It 
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to the ttN NKK reaction is given by the diagram 


Then all reactions are related to four isospin combina¬ 
tions parameterized in as: 


N 


K 

K (B8) 

(T(7r“ A'*"''' - 

K+TP) = 

N 

tT(7r°A° - 

4 A+E-) = 

be parameterized as 

AP - 

A+E°) = 


a{^aNb ^ N,KdK,) = cl.l2l 


■l_f£- 

1.86 

■fo' 

s . 


- s . 


,(B9) 


where = rriN 2mK = 1.929 GeV and c is a channel 
dependent isospin coefficient. The latter ones are col¬ 
lected Table m In Figs. CZI and na we confront these 
parameterizations with the available experimental data. 


0.004959(Vs- 1.688)°™® 
(v^- 1.725)2 -^0.008147 ’ 
0.006964(Vs- 1.688 )°®i 4° 
(v^- 1.725)2 -^0.007713 ’ 
0.002053(Vs- 1.688)°°®®° 
{y/s- 1.725)2 -^0.005414 ’ 
0.3179(Vs- 1.688)°°°2® 
(Ts- 2.675)2 -^44.88 ’ 

^ ' (yi- 1.725)2-t 0.003777’ 

For other reactions one finds 


2. Reactions of tt -|- A 


The cross sections for ttA ^ reactions we take 

from Ref. 


cr(7r"A++ ^ K+A) 


0.00983(Vs- 1.613)°™° 
(^s- 1.720)2 -^ 0.004852 ' 
(BIO) 


The other channels leading to A production are obtained 
from isospin symmetry 

(j{'K~A) = cr(7r''"A~ ^ iF°A), 

(T(7r°A+ ^ AT+A) = (T(7r°A° ^ iF°A) 

= ^a(7r-A++ ^ K+A), 
cr(7r+A° ^ AT+A) = cj{'k~ Ar°A) 

= ia(7r-A++ ^ K+A), 

(Bll) 


(T(7r°A°’''' - 

A+E+) 

= (T(7r°A“ — 

> A°E-) = 0, 

cr(7r“'"A+ — 

^ A+E+) 

= (T(7r~A° — 

> A°E-) =0, 

cr(7r+A“ - 

.4 A°E°) 

= a{'K~ 

^ A+E°), 

cr(7r°A+ — 

^ A+E°) 

II 

O 

> 

o 

i 

A°E°) 




^ A+E-), 

CT(7r“A+ - 

A°E°) 

= CT(7r^A° — 

^ A+E°), 

CT(7r“A+ — 

> A+E-) 

= CT(7r^A° — 

^ A°E+) 



= |a(7r-A++ ^ A+E°), 

cr(7r°A+ — 

^ A°E+) 

= (T(7r°A° ^ 

A+E-), 

<7{'K~ - 

A°E+) 

= Cj{'K^ A~ - 

. A+E-), 


(B12) 

We also include reactions ttA ^ NKK. These are 
taken with the cross section parameterized as in (IB9II . 
The corresponding isospin constant c follows from the 
analysis of the diagram (IB8II . where the incoming nu¬ 
cleon line is replaced by the A, and depends on incoming 
and outgoing isospin via Clebsch-Gordan coefficients (cf. 
1^1. The coefficients are collected in Table EH 


For ttA ^ EAT reactions one assumes in that 
the cross section has a contribution from the isospin 1/2 
channel only. The / = 3/2 contribution was found small. 


TABLE V: Isospin coefficient in the parameterization 1B9II of 
hN —> NKK reactions 


reaction 

c 

reaction c 

TT+P —> pA+A° 

1 

Tv~p —> pK^K~ 1 

TT+n ^ pA°A° 

2 

TT-p^nK°K° 2 

TT^n —> pK^K~ 

2 

■K~p nK^K~ 2 

Ti'^n nK^K^ 

1 

TT°n^pK°K- 2 

Tv°p pK°K° 

1 

2 

TT°n^nK°K° f 

7r°p —> pK^ K~ 

1 

2 

7r°n —> nK^ K~ i 

n°p nK+K° 

2 

Ti-n^nK°K- 1 


3. Reactions of A -|- A 

Gross sections for reactions of the type AA —> NAK 
and AA —> AEA are parameterized in |^. For A A ^ 
NAK reactions we have 

cr{pp pAK^) = a{nn —> nAA°) 

= 1.879 l) 

\SkNA ) \ S J 

a{pn nAK^) = a{np pAA°) 

= 2.812 l) (^^)^®%13) 

\SkNA j \ S j 

where skna = (wa -I- toat -I- rn^)^ . These parameter¬ 
izations are confronted with the available experimental 
data in Fig. ^1 For AA ^ ASA reactions we use 
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FIG. 17: Total cross section of the —> pK^K^, n p ^ nK°K^, tt p ^ nK^K , and n p ^ pK^K reactions: solid line 

calculated with eq. dots are the experimental data from Ref. 




FIG. 18: Total cross section of the rr+n —> pK^K~, rr+n —> 
pK°K°, and rr+n ^ nK^ reactions: solid line calculated 
with eq. llB^ . dots are the experimental data from Ref. (^. 
Note that isospin symmetry of the diagram IIB8II predicts that 
cross sections for Tv'^n pK^K~ and the one for 7r"''n —> 
pK°K° multiplied by 4 are same. 


TABLE VI: Isospin coefficient in the parameterization 1IB9II 
of ttA —> NKK reactions 

reaction c reaction c 

7r+A+ ^ pK+K° 2 


7r° A++ 

^ pK+K° 

3 

7r+A° 

^ pK°K° 

1 

7r+A° 

pK+K- 

1 

7r+A° 

^ nK+K° 

2 

7r° A+ 

^ pK°K° 

1 

7r° A+ 

pK+K- 

3 

7r° A+ 

^ nK+K° 

1 

7r"A++ 

pK°K° 

3 

7r"A++ 

pK+K- 

3 

7r“A++ 

nK+K° 

0 


tt+A- ^ nK+K- 3 
tt+A" ^ nK°K° 3 
tt+A" ^ pK°K- 0 
7r° A° ^ pK°K- 1 
Tv° A° ^ nK°K° 3 
■K° A° ^ nK+K- 1 
7r-A+ ^ nK+K- 1 
7r-A+ ^ nK°K° 1 
7r“A+ ^ pK°K- 2 
7r° A" ^ nK°K- 3 
7r-A° ^nK°K- 2 


FIG. 19: Total cross section of the pp pAK+, pn pAK^ 
reactions: solid line calculated with eq. (ini3, dots are the 
experimental data from Ref. 


(j{pp pTPK+) = a{nn -i- uTPk'^) 

2.753 


= 5.321 


- 1 


^SkNS 


(^) 


8.510 


a{pp ^ nT,+K+) = a{nn—> pT, K^) 

2.743 


= 1.466 


- 1 


-K°) 

/ Sif ATE 


.SifATE 

a{pn p'E,-K+) = a{np nY]+K^) 

2.782 


= 11.02 


- 1 


\ s J 

:+K°) 

/ Skate \ 7.674 




.Skate 

a{pn nT,^K+) = a{np pTPK^) 

2.773 




= 6.310 


- 1 


/ Skate \ 7-820 




^sknt, 

a{nn vAl-K+) = a{pp p'S+K^) 


s 

.'0\ 




= 7.079 


- 1 


Skate 


2.760 


(^) 


8.164 


(B14) 


where Skate = (we + toat + ttikY■ These parame- 
terizations are compared with the experimental data in 
Figs. 1201 ini HOI 

Next we proceed to reactions of the type NN 
NNKK. Following we will use the same isospin- 
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FIG. 20: Total cross section of the pp —> pTi'^K^, pp —> 
nTj'^K'^ reactions: solid line calculated with eq. IIB14L dots 
are the experimental data from Ref. 


TABLE VII: The isospin coefficients in parameterizations of 
the reactions NN AAK and NN AEA i5l7l5l8lBl9ll 



reaction 

CA 


reaction 

C-E 


reaction 

CS 

PP 

^ A+AA+ 

1 

3 

pp - 

A+ S°A+ 

1 

6 

pn - 

A+ eO K° 

1 

6 

pn 

^ A+AAO 

1 

3 

pp - 

A+ S+A“ 

1 

3 

pn - 

-* A° T.+K° 

1 

3 

pn 

^ AO AK° 

1 

3 

pp - 

A++sOaO 

1 

2 

pn - 

A++T.-K° 

0 

nn 

AO AK° 

1 

3 

pp - 

-^AO S+A+ 

0 

nn • 

^ A° S-A+ 

1 

3 

nn 

A-AK+ 

1 

pn ■ 

AO j:°k+ 

1 

6 

nn ■ 

A- E°K+ 

1 

2 




pn • 

A+ S-K+ 

1 

3 

nn • 

-^AO eO^o 

1 

6 




pn ■ 

A- T.+K+ 

0 

nn • 

-^A- E+K° 

1 







nn ■ 

A+ E-aO 

0 



FIG. 21: Total cross section of the pp —> pTPK^, pn 
pT,~ reactions: solid line calcnlated with eq. IIB14II . dots 
are the experimental data from Ref. 


of reactions pp and pp K'^ are 

parameterized as follows 


(j{pp A ) = 6.166 

a{pp A++E-A:+) = 10.00 


s 

SAAK 

S 

SAEiC 


1 2.842 

Saak' 


S - 


(B17) 


1 2.874 

sask' 


s 


(B18) 


where saak = ijnA + rriA + rriK)^ and sat.k = (wa + 
mg + itik)^ ■ The other reactions channels are related to 
these two by isospin coefficients ca and ce 


averaged cross sections 

a{NN NNKK) = 1.5 (^1 - — j j (B15) 


= (2mAr+ 2mA')^ 

for all possible isospin channels 

pp- 

-> ppK^K 

pn - 

ppK^K~ 


pp- 

ppK°K° 

pn - 

nnK'^K^ 


pp- 

pnK+K° 

nn - 

nnK^K~ 

(B16) 

pn - 

-> pnK^K~ 

nn - 

-4 nnK^K^ 


pn - 

-4 pnK°K° 

nn - 

pnK^K~ 



Cross sections for processes NN ^ AAK and NN ^ 
AT,K were calculated in Ref. [^. The cross sections 



a{NaNb AcAKd) = CAcr{pp A++AAr°) 

a{NaNb AJ^dKe) = CY,(T{pp A++E“Ar+) 

(B19) 

This coefficients are collected in Table ivrn 


4. Reactions of A + A 


Reactions AN NYK {Y = A, S) are parameterized 

in 1^ in terms of two cross sections 


cr(A++n ^ pAK+) = 8.337 
a{A-p nY-K+) = 52.72 


- 1 


n 2.227 


/ 

[™] 


-1 


1 2.799 


[™] 


2.511 

(B20) 

6.303 

(B21) 


The other reaction channels are related to these two by 
the isospin coefficients ca and ce 

a{AaNb NcAKd)=CA cr(A++n ^ pAK^) 

a(AaNb NcYdKe) =ce a{A~p nE“Rr+) 


FIG. 22: Total cross section of the pn —> pTPK° reaction: 
solid line calculated with eq. dots are the experimental 

data from Ref. l3^ 


(B22) 
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TABLE VIII: The isospin coefRcients in parameterizations of 
the reactions AN —> NYK IIB20IB21IB22II 


For reactions of the type AN NNKK, all channels 


^AA — 


reaction 


CA 

reaction 


ce 

reaction 

Ce 

^AA — 

A++n — 

^ pAK+ 

1 

A++P ^ 

- pS+A+ 

0 

A+n - 

-^pT.-R+ 

1 

3 

A+ p ^ 

pAK+ 

1 

3 

A+ p ^ 

p’e+ro 

1 

3 

AOp- 

■* pS-A+ 

1 

3 


A+ n ^ 

pAK° 

1 

3 

A++n — 

>pE+A0 

1 

A+n - 

nS0A+ 

1 

6 


A° 

pAK° 

1 

3 

A+ p ^ 

pS0A+ 

1 

6 

AOp- 

-+ nT,0R+ 

1 

6 

AO) _ 
^AA — 

A+ n 

nAK+ 

1 

3 

A++n — 

>pS°A+ 

1 

2 

A°n- 

nT.0R0 

1 

6 

A° 

nAK+ 

1 

A^" p —>■ 

nS+A+ 

0 

A~p - 

-> uSOrO 

1 

= 

A~ p —> 

nAK° 

1 

A++n — 

> nS+A+ 

0 

A°rt- 

■^pT.-R0 

0 

^AA — 

A° 

nAK° 

1 

3 

n 

ptPro 

1 

6 

A~p - 

■*pT.-R0 

0 




AO p^ 

pT.0r0 

1 

6 

AOn- 

nY,-R+ 

1 

3 





n 

nT,+R0 

1 

3 

A~p - 

-+ nY,-R+ 

1 





AO p^ 

nS+AO 

1 

3 

A~n - 

nH-RO 

0 

^AA — 


(9) 

^AA 


cr(A++p ^ a++e°a:+) 

n 2.748 


= 10.30 


- 1 


SASK 

cr(A++n ^ A++E"A:+) 


[^] 


9.321 


= 10.33 


1 2.743 


- 1 


fSASif 18-915 


SAS/f 

cr(A+p ^ a+e°a:+) 

1 2.759 


= 10.62 


- 1 


SASK 

a{A+p A°E+A:+) 

n 2.830 




= 0.647 


- 1 


_ SAS/f 

a{A°p^ A+T,-K+) 

n 2.843 




= 2.128 


- 1 


SAS/f 

cr(A+n ^ A+E-A:+) 

n 2.757 


[^] 


= 10.57 


SASK 


- 1 


[^] 


10.20 


3.862 


5.986 


10.11 


A++P - 

-4 ppK+KO 

A^p - 

-> puK^Ro 

A++n - 

ppK^RO 

A°p - 

nnK^RO 

A++n - 

-> ppK'^K~ 

A°p - 

PPK0k~ 

A++n - 

pnK'^RO 

A°p - 

■> pnK~^K~ 

A+ p - 

-> ppKORo 

AOn- 

nnK^Ro 

A+ p - 

-> ppK^K~ 

A^n- 

pnK0K~ 

A+ p- 

-t pnK^RO 

A^n- 

-> nnK'^K~ 

A+ n - 

priK^Ro 

A~p - 

nnKORo 

A+ n - 

nnK^RO 

A~p - 

pnK0K~ 

A+ n - 

PPK0K~ 

A~p - 

nnK^K~ 

A+ n - 

pnK^K~ 

A~n - 

nnK0K~ 


(B23) 


will be taken with the averaged cross section ior NN ^ 
NNKK [^. 

Next, we consider the family of processes AN AYK 
(Y = A, E). In Ref. the following ten reactions were 
parameterized: 


^AA 


cr(A++p ^ a++ar:+) 

n 2.303 


= 2.704 


- 1 


( 2 ) 

^AA 


SAAK 

a{A+p A+AK+) 

2.350 




= 2.917 


- 1 




5.551 


6.557 


'I’aa = (7(A'^n ^ A^AK'^) 


= 0.312 


- 1 


2.110 


[“]''i®5, (B24) 


^(10) 

lAA 


= cr(A+n ^ A°E°R:+) 

n 2.846 


= 1.112 


- 1 


5.943 ^225) 


The other isospin channels can be expressed in terms of 
these ten cross sections as it is given in Table im 


5. Reactions A + A 

Reactions A A ^ AYK have small cross sections and 
thus contribute marginally to K'^ production. Neverthe¬ 
less, they are included in the calculation with parameter¬ 
izations from 1^. 


^AA = o-(A+A++ ^ a++ar:+) 


AA 

= 1.154 

AA = ^(A°A 
= 0.881 

=.(3) 


- 1 


2.149 


/SAAif 




SAAK ) \ S ) 

oa++ . a+ar:+) 

S 7-977 


U 


AA 


,saak 

= ct(a°a+^ a°ar:+) 


Vs/ 


= 0.291 


- 1 


SAAK 


2.148 


^ SAAA- ^ 7.934 


(B26) 


A4) 

'AA 


(J\\ = cr(A++A° ^ A++E"R:+) 

^ 2.953 

= 3.532 ( —^-1' 


SAT.K 


^SAE^y2.06 


= cr(A++A° ^ a+e°r:+) 


= 2.931 


- 1 


5ASiC 


2.952 


(^) 


12.03 
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TABLE IX: Cross sections of reactions AN —> AYK (Y = A, E) in terms of cross sections 1IB24IB25II 

AN AAA 


reaction 


a reaction 


(T reaction 


a reaction 


(T reaction 


A++P ^ A++AA+ 
A+ n ^ A+ AK° 
AO A- \K+ 


<^AA ” 

4^(3) A- 
s'^AA ^ P 


*A++KKO 
> AO AK+ 
A° AKO 


'^AA 

4^(3) 


A++n ^ A+AA+ 
AO A+ARO 
A- A-AK+ 


|4a 
4i AV- 

4a A~n- 


> A++AA° A+p ■ 

. A“ AK+ AO n 

. A- ARO gW 


. A+AA+ 
.AOAAO 


AA ■ 


ASA 


reaction 


reaction 


reaction 


a reaction 


a reaction 


A++P 


a++e”a:+ 

^AA 

A+p 


A++E“if° 

3^(10) 
4 ^AA 

A+n ^ A°'S°K+ 

^(10) 

°AA 

A°p 


A“E+if+ 

4 ^AA 

A~p —»■ 

A+S-rO 

0 

A++P 


A++E+A'“ 

^AA 

A^p 


A+'S°K+ 

•^AA 

A+n ^ A++S-if“ 

3^(7) 

4 "^AA 

A°n 


a°yPk° 

'Taa 

A~p —»■ 

A-Y.+ R0 

^AA 

A++P 


A+S+A'+ 

3^(7) 

4 ^AA 

A+p 


A+S+K° 

•^AA 

A+n ^ A~S+if+ 

0 

A“n 


A“S“iC+ 

•^AA 

A^p 

A~Y.°R+ 

^AA 

A++n 


a++s“a“ 

^AA 



A++'S-K+ 

4 ^AA 

A°p A+SOrO 

^(10) 

^AA 

A^n 


A+E-iC“ 

•^AA 

A~n —»■ 

A~T.°R° 

•^AA 

A++n 


A+E“A'+ 

3^(10) 

4 "^AA 

A“*"p 


A°^+K+ 

•^AA 

A°p A + E-iC+ 

‘^AA 

A^n 


A-Y.+K° 

4 ^AA 

A~n 

A-T.-R+ 

•^AA 

A++n 


A+S+if® 

4 ^AA 

A+n 


A+S°K° 

^AA 

A°p A“E+if“ 

'^AA 

A“n 


A^S“if+ 

3 _(10) 

4 ^AA 

A~n 

A°S-R° 

3^(7) 

4 '^AA 

A++n 


A++E-A'+ 

•^AA 

A+n 


A+E-K+ 

•^AA 

A“p ^ A°E“if+ 

‘^AA 

A~ p 


A°E“if° 

3^(10) 

4 ^AA 




A++n 


a“s+a:+ 

0 

A'*'n 


A°S+if“ 

^AA 

A°p A+ + S-K° 

0 

A~ p 


A“E-if+ 

4 ^AA 





^(6) 

^AA 


A7) 

’AA 


cr(A'A+ ^ a°e-a:+) 


5.861 


(— 

\SA-SK 



S J 


a(A-A° A-S-K+) 

7.047 (— -l) (B27) 

\SA^K / Vs/ 


The other isospin channels can be expressed in term of 
these seven cross sections as it is given in Table El 
Another group of reactions are A A —s- NNKK. We 
use the same parameterization (IB15II for all isospin chan¬ 
nels 


A++A+- 

ppK^K^ 

A-A++ - 

-4 pnK^K'^ 

A++A0- 

-4 ppK^K^ 

A-A++ - 

-4 nnK'^K^ 

A++A0 - 

-4 ppK~^K~ 

A-A++ - 

-4 ppK^K~ 

A++A° - 

-4 pnK'^K^ 

A-A++ - 

-4 pnK~^K~ 

A+ A+ - 

-4 ppRORO 

AO A° - 

4 nnK^K^ 

A+ A+ - 

-4 ppK^K~ 

AO AO - 

4 pnK^K~ 

A+ A+ - 

-4 pnK^K^ 

AO AO - 

4 nnK~^K~ 

A+ A° - 

4 pnK^K^ 

A-A+ - 

-4 nnK^K^ 

A+ AO - 

4 nnK~^K^ 

A-A+ - 

-4 pnK^K~ 

A+ AO - 

4 ppK^K~ 

A-A+ - 

-4 nnK^K~ 

A+ AO - 

4 pnK~^K~ 

A-AO - 

4 nnK^K~ 


6. Reactions of tt -|- tt 

The source of information about tttt KK reactions 
is the inelastic tttt scattering which was intensively stud¬ 
ied experimentally ^3- ^ characteristic feature of the 
TTTT scattering is very rapid rise of the phase shift and the 
inelasticity parameter at energies about 1 GeV. Such a 
strong energy dependence is caused by the /o resonance 
with a pole around 980 — il4 MeV (cf. Ref. ^3)- The 



FIG. 23: The cross section of the 7r"'"7r ^ K'^K reaction: 

solid line represents parameterization dashed line is 

calculated as in Refs. [3^. l44). The experimental points are 
from 


inelasticity parameter translated into the cross section 
of the 7r+7r“ —> K~^K~ reaction is depicted in 

Fig. US We checked, that the velocity-averaged cross sec¬ 
tion is reproduced almost exactly if we used a constant 
for a TTTT —!■ KK 

CT(7r^7r“ —!■ K~) = CT(7r“7r“'' ^ = 1.6 mb , 

(B29) 

shown in Fig. [221 by the solid line. ,i,From isospin sym¬ 
metry, the other tttt channels have cross sections 

cr(7r°7r° ^ K+R-) = cr(7r°7r° ^ K°K°) (B30) 

= ^ a{TT^TT~ K~) = 0.64 mb, 

5 

^ K+K°) = -> K°K-) (B31) 

= — trlTT^Tr” ^ K'^K~) = 1.92 mb . 

5 

Note that this implementation is different from the tttt 
cross sections used in The latter is shown in 

Fig. 1221 by the dashed line. 
























20 


TABLE X: Cross sections of reactions A A ^ AY K (Y = A, E) in terms of cross sections IIB26IB27t 


AA - 


AAK 


reaction 


a reaction 


a reaction 


a reaction 


O' reaction 


A+ + A+ ^ A++AK+ 

A++A° ^ A++AK° 

0 

A++A° ^ A+AK+ 


A+ A+ ^ A++AiC“ 

1 =(1) 

3 '^AA 

A+A+ ^ A+AK+ 

9 °AA 

A+ A“ ^ A+ AK° 

A+ A“ ^ A“ AK+ 

^(3) 

A+ + A“ ^ A+AK° 

0 

A++A“ ^ A“ AK+ 

0 

A“ A” ^ A°AK° 

9 °^AA 

A° A“ ^ A" AK+ iff// 

A- A+ ^ A° AK° 

•^AA 

A- A+ ^ A-AK+ 

0 

A- A° A~ AK° 

'^AA 




AA ^ AYK 


reaction 


0 

reaction 

cr 

reaction 



a 

reaction 

0 

A++A+ 

^ A++S°A'+ 

2'^AA 

A+A+ ^ A+E+E'° 

g'^AA 

A"A++ 


A+E-A'+ 

'^AA 

A-A+ ^ A°E°E'° 

^AA 

A++A+ 

^ A++Y+K° 


A+A+ ^ A++E-E'+ 

3^AA 

A-A++ 


A°S+E''’ 

'^AA 

A-A+ ^ A°E-E'+ 

•^AA 

A++A+ 

A+T,+K+ 

0 

A+A+ ^ A°S+E'+ 

0 

A“A++ 


A°S°E'+ 

^AA 

A-A+ ^ A+S-E'° 

0 

A++A° 

A++T,°K° 

^ AA 

A+A° ^A+S“7L° 

6^AA 

A"A++ 


A++E-E'‘’ 

0 

A-A+ ^ A-E+E'° 

•^AA 

A++A° 

A+T,°K+ 

'^AA 

A+A° ^ A+E-E'+ 

2^(5) 
3 " AA 

A~A++ 


A-E+1L+ 

0 

A-A+ ^ A-E°E'+ 

'^AA 

A++A° 

A+'£+K° 

'^AA 

A+A° ^ A°S+7L° 

2^(5) 

3"AA 

A° A° 


A°T,°K° 

IS^AA 

A-A° ^ A-TPK° 

2‘^AA 

A++A° 

A++Y-K+ 

^ AA 

A+A° ^ A°S°E'+ 

6^AA 

A° A° 


A°T,-K+ 

9*^AA 

A-A° ^ A-E-E'+ 

^AA 

A++A° 

A°'£+K+ 

0 

A+A° ^ A++E-E''’ 

0 

A'’ A° 


A+T,-K° 

0 

A-A“ ^ A'’E-E''’ 

0 

A+ A+ 

A++T,°K° 

e'^AA 

A+A° ^ A-E+E'+ 

0 

a“ A° 


A-T,+K° 

3^AA 

A++A++ ^ A++S+if+ 

0 

A+ A+ 

A+T,°K+ 

18 "^A A 

A-A++ ^ A+S°E'‘’ 

°^AA 

A“ A° 


A-T,°K+ 

fi*^AA 

A- A- ^ A-T,-K° 

0 


7. Reactions of p + p 


The cross section for these reactions is not known, be¬ 
cause there are no data to which calculations could be 
compared. It has been calculated in as a kaon ex¬ 
change process. We adopt this result here and parame¬ 
terize the cross section as 


cr(p+p 


K+K-) = (T{p+p- K°K°) (B32) 


^ Pcm 
327r p^m 


exp 


a/s - 2mp \ 
IGeV ) 


29.11 mb, 


den. The reaction is dominated by / resonance 

3 


a{TTp (j) —> KK) — 


At: s ^ 


(^) (^) 


Pi (s-+ sr2ot(\/s) 

(B35) 

where pi is the momentum of incoming particles in the 
centre-of-mass system and 


qo = PcmiM^, niK, rriK) , q = Pcm(s, mx, rriK) , 

ko = Pcm(AL|, 771^, TOp) , k = Pcm(s, 771^, 777p) , 

and the total width is determined from the partial widths 


where _p^(s), and Pcmi^) 8’'^® momenta of kaons and 
rhos in the centre of mass frame, respectively. The expo¬ 
nential form is chosen such that it reproduces very well 
the matrix element calculated in [T^. 

Other pp channels are determined from isospin sym¬ 
metry, in analogy to tttt scattering 


0 

0 

-ff K"*K-) = cr(pV° ^ K^K^) 



= ^cr(p+p“ ^ K+K-), 

(B33) 

a(p+p° - 

K+K°) = a{p-p° K°K-) 



= ^o-(p+p" ^ K+K-). 

(B34) 


8. Reactions of n + p 

Production of a ATpair in this reaction is suppressed 
because the s-wave component of cross section is forbid- 


0 

rtot('\/s) = + r 

Kq 


- ) 

-3 

qo 


(B36) 


The partial widths for decays of / into K'^K~ and K^Kg 
are different, but since we keep isospin symmetry in all 
other reactions throughout this work, we shall use 




(B37) 


for all / = 0 channels. 

For the reaction of Trp into KK*, the leading or¬ 
der diagram in chiral counting is the contact Weinberg- 
Tomozawa term |4^ . We calculated the cross section 


(T 


(B 


o-x 


9 Pcni{s,mx,iTlK*) | ,^^(J,0)|2 j _ q 
647rs Pcm(s,m.^,mp) 

^ Pcmis,mK,mK*) ^ /^(0,0) rjVf(L0)j^ 
647rs pcni{s,m^,mp) V ^ ’ 


using unitarized amplitudes from where I 

and S are isospin and strangeness quantum numbers in¬ 
dicating a particular reaction channel. The cross sections 
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FIG. 24: Cross sections for reactions up —> KK* with given 
isospin and the mixed terms ax 


are shown in Fig. We realized that, to a very good 
approximation, only the J = 1 component of the cross 
section effectively contributes and can be replaced by a 
constant. From combining the involved isospin states we 
obtain relations 


+ 

4 

1 

K+K*-) = (j{p--K+ 

K°K*°) = 0.2 mb, 

1 

q 

+ 

K'^K*~) = (j{p'^TT~ 

^ K°K*°) = 0.2 mb, 

O 

o 

1 

^ K+K*-) =cr(p°7r° ^ 

, x°K*°) = 0 

1 

o 

+ 

-4 K+K*°) = a{p\° ^ 

K°K*°) = 0.4 mb. 

cr(p°7r^ - 

-4 K+K*°) = a{p°TT- - 

, K^K*-) = 0.4mb. 


(B38) 


10. Reactions of F + tt 


These reactions produce kaon and a cascade. The cross 
sections are unknown, so we just use a parameterization 


ncm(s,TOH,TOic) 
(7 G , , 

Pcm{s,mY,m^) 


10 mb, 


where the constant C results from the Clebsch-Gordon 
coefficients of isospin adding. In the following we just list 
the parameterizations and the corresponding constants 
C: 


A7r+ 

o 

[I] 

-H 

T 

Att — 


1 (B42) 

2 

AttO 

^ A+S- , 

AttO ^ 


E+ttO ■ 

A+SO , 

E-TfO 

-^K°E- 

1 

^ 2 

E“''7r“ 

^ A+S- , 

E-7r+ 


5 

12 

E0 7r+ 

o 

[I] 

+ 

T 

E^tt- 

-^K°E- 

^ i (B43) 

o 

o 

A+S- , 

EOttO - 


1 

^ 3 

E Tr’*' 

^ A+S- , 

E+7r“ 

^K°E^ 

5 

12 

11. 

Total cross 

sections of KN reactions 


These parameterizations are taken from |48j. 


{ 12.4 

1.09 + 14.5p_R' 
18.64-0.5pif 


PK < 0.78 

0.78 <PK < 1.17 (B44) 
1.17 <PK < 2.92 


9. Decays of K* 


Kaons can be also produced in decays of K*. The 
decay rate is simply given by Tk*Pk* where F^f. = 
50.8 MeV is the width and pK* the density. Individ¬ 
ual channels must be multiplied by appropriate Clebsch- 
Gordan coefficients. These are 


K*+ K+Tr° , K*° K°Tr° 

K*+ K°tt+ , K*° K+tt- 



(B39) 


The inverse reaction destroys kaons. The cross section is 
derived as ^3 


{ 15.5 

4.19 -I- 14.5p_R- 
22.78 - 2ApK 


Pk < 0.78 

0.78 <PK < 1-1 (B45) 
11 <Pk < 2.3 


where pk is kaon momentum in the lab frame. In Fig. EH 
we illustrate the quality of the paremeterization 


APPENDIX C: FEED-DOWN FROM 
RESONANCE DECAYS 

We list in Table IXII the average numbers of pions 
which are obtained by decays of individual resonances. 
Resonances which feed-down into kaon, A and E produc¬ 
tion are listed in Table mu 


t{Ktt K*) = 


TrFic. Ak* 






,mK) 


(B40) APPENDIX D: STRANGENESS CONTENT OF 

THE INITIAL STATE 


where 


1. Strangeness mnltiplicity in NN collisions 


Ak* (Vs) = 




(s - 


(B41) Strange particles are also produced in the initial 
nucleon-nucleon interactions. We shall assume that their 
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FIG. 25: The total K^p cross section as given by IIB44t in 
comparison with the experimental data from 


TABLE XI: List of resonances decaying into pions together 
with the average numbers of pion species produced by their 
decays. 




)r 


{tt+)r 

(tt )r 

V 

0.28 

0.28 

iVO(1440) 

0.175 

0.725 


1 

0 

iV+(1440) 

0.175 

0.725 


1 

1 

AO (1440) 

0.725 

0.175 

p~ 

0 

1 

K*+ 

2/3 

0 

OJ 

0.91 

0.91 

A*o 

0 

2/3 

A++ 

1 

0 

K- 

0 

2/3 

A+ 

1/3 

0 

A*o 

2/3 

0 

AO 

0 

1/3 

A(1405) 

1/3 

1/3 

A- 

0 

1 

A(1520) 

0.217 

0.217 

A++ 

0 

1 

E+(1385) 

0.94 

0 

A+ 

0 

1/3 

E0(1385) 

0.06 

0.06 

AO 

1/3 

0 

E-(1385) 

0 

0.94 

A- 

1 

0 


0 

1 

^+(1440) 

0.725 

0.175 

n 

0 

0.322 


TABLE XII: List of resonances decaying into kaons, lambdas 
or sigmas together with average numbers of decay products. 



{K+)r {K-)r {A)r (E+)fl (E-)k 

K*+ 

1/3 

K*o 

2/3 

K*- 

1/3 

K*° 

2/3 

eO 

1 

A(1405) 

1/3 1/3 1/3 

A(1520) 

0.11 0.143 0.143 

E+(1385) 

0.94 0.06 

e0(1385) 

0.88 0.06 0.06 

E-(1385) 

0.94 0.06 


1 


1 

n 

0.678 1 


TABLE XIII: Multiplicities of negative hadrons i n pp, pn , and 
nn collisions interpolated from data compiled in 14911^ . The 
averaged values for nucleon-nucleon collisions are calculated 
according to eq. ins. 



{h )pp 

')pn 

)nn 

{h )nn 

Au-I-Au @ 11.6 AGeV 

0.697 

0.900 

1.390 

1.044 

Pb-bPb @ 30 AGeV 

1.265 

1.728 

1.965 

1.741 

Pb-bPb @ 40 AGeV 

1.482 

1.899 

2.182 

1.937 

Pb-fPb @ 80 AGeV 

2.025 

2.434 

2.725 

2.477 

Pb-tPb @ 158 AGeV 

2.611 

3.058 

3.311 

3.081 


initial multiplicities are given by production of kaons in 
NN reactions in vacuum. Since in our approach kaon 
production is calculated explicitly and species with S' < 0 
are populated statistically, we will only need input on pri¬ 
mordial production. Data on kaon multiplicities in 
pp collisions were summarized in . In order to obtain 
the average multiplicity in nucleon-nucleon collisions, we 
also need multiplicities in nn and pn collisions. It was 
observed that in TViV collisions multiplicity of charged 
kaons does not depend on isospin of the projectile/target 
and we can represent {K'^)nn with {K^)pp [^1^ 

{K+)nn = {K+)pp = {K+)nn = {K+)pn . (Dl) 


2. Density from multiplicity 

What we actually need for the simulation is the den¬ 
sity, not the multiplicity. We will estimate the former by 
equating 


^ (D2) 

Ph- {h }nn 

Data on multiplicity of negatively charged hadrons from 
NN interactions were collected by Gazdzicki and Rohrich 
in 1^. They collected h~ yields from pp and pn colli¬ 
sions at different energies and estimated them for nn col¬ 
lisions. The primordial yield (per nucleon-nucleon pair) 
for nucleus-nucleus collisions is then obtained as 

2 

{h )pp 

+ 2| (l - I) {h-)pr. + ■ I)' 

This relation is formulated for nuclei with atomic number 
Z which contain A nucleons. 

The values for h~ yields are not measured at all en¬ 
ergies we need. Thus we use data collected in mil 
and linearly interpolate as a function of ^/s. The com¬ 
pilation is displayed in Table IXITII The resulting ratios 
{K'^)nn/ {h~) NN are listed in Table nn 
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3. Composition of multiplicity data and feed-down 

The major part of the h~ multiplicity are pious, 
a smaller fraction comes from negative kaons. Very 
small contribution is due to antiprotons. There is feed- 
down from resonance decays mainly into pion multiplic¬ 
ity. Feed-down to antiprotons is neglected because an- 
tibaryons are populated scarcely in general at these en¬ 
ergies. Feed-down to K~ and from K* is taken into 
account. 

This is implemented by using the effective densities of 
pions, kaons and antiprotons when determining the h~ 
density at the left-hand-side of eq. iTH^ 

Ph~ = P#- + Pk- + Pp ■ (D4) 

The effective densities, denoted by tildes over subscripts, 
include the actual densities and pions (kaons) which can 
be produced by decays of all resonances in the system. 
These additions are calculated according to the same pre¬ 
scription as the feed-down in final state, see Section ITTfI 
and Appendix 10 By using the same strategy we keep 
the simulation internally consistent. 

4. Algorithm to obtain the initial state 

Technically, the initial state is obtained iteratively: 

1. Set densities of all strange particles to zero. Set 
temperature and chemical potential to some “rea¬ 
sonable” value. 

2. Determine energy density and I 3 density contribu¬ 
tions stored in kaons, and psic, respectively. 
This is done by first determining phase space occu¬ 
pancies via 



where j = K'^, . Then 

X (7x+ exp(/r3/2T) -k7A:o exp(-^3/2T)) , (D6) 


and 

P3K =-^{Pk+- Pro) ■ (D7) 

Subtract them from the total e and ps, i.e., calcu¬ 
late 

e = e — ek (D8) 

P3 = P3 — P3K (D9) 

3. Calculate ps = Pk+ + Pk° ■ 

4. /^From the new e, pB, ps, ps calculate temperature, 
chemical potentials and the suppression factor 75 . 

5. From T, pB, P 3 , is calculate p*- -I- -I- pp. 

This also includes resonance feed-down to pions 
and kaons according to the given temperature and 
chemical potentials. 

6 . Get A'+ density as 

{K+)nn , , 

P^+ = 77 —t- (P*- PB, P3, is) 

[n /jvAT 

+ Pp{T, Pb, P 3 ,7s) + Pk~ PB, P 3 ,7s)) (DIO) 

7. Determine 

Pro = Pr+ exp(-p 3 /T) (Dll) 


8. Calculate again ps = pk+ + Pro- 

9. If the obtained ps differs from that determined in 
step 13 go to step 13 If they agree, the routine has 
converged. 
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